Integrability of the Minimal Strain Equations for the Lapse and Shift in 3 + 1 

Numerical Relativity 
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Brady, Creighton and Thorne have argued that, in numerical relativity simulations of the inspiral 
of binary black holes, if one uses lapse and shift functions satisfying the "minimal strain equations" 
(MSE), then the coordinates might be kept co- rotating, the metric components would then evolve on 
the very slow inspiral timescale, and the computational demands would thus be far smaller than for 
more conventional slicing choices. In this paper, we derive simple, testable criteria for the MSE to be 
strongly elliptic, thereby guaranteeing the existence and uniqueness of the solution to the Dirichlet 
boundary value problem. We show that these criteria are satisfied in a test-bed metric for inspiraling 
binaries, and we argue that they should be satisfied quite generally for inspiraling binaries. If the 
local existence and uniqueness that we have proved holds globally, then, for appropriate boundary 
values, the solution of the MSE exhibited by Brady et. al. (which tracks the inspiral and keeps the 
metric evolving slowly) will be the unique solution and thus should be reproduced by (sufficiently 
accurate and stable) numerical integrations. 

PACS numbers: 04.25.Dm, 04.30.Db, 04.70.-s 



I. INTRODUCTION 

Brady, Creighton and Thorne 0| have shown that 
post-Newtonian techniques are unable to treat the late 
stages of inspiral of black- hole binaries with sufficient ac- 
curacy to satisfy the needs of gravitational-wave data 
analysis. The only viable way to perform these com- 
putations appears to be numerical relativity. The chal- 
lenge of performing such computations is called the In- 
termediate Binary Black Hole (IBBH) problem: the 
late inspiral is intermediate between the early inspiral 
(where post-Newtonian techniques are successful) and 
the merger (where numerical relativity must be used). 
No researchers have yet tackled the IBBH problem by 
numerical relativity, but it is a high-priority, near-future 
challenge [0. 

If conventional numerical relativity techniques are used 
for the IBBH problem, then the coordinate system will 
be asymptotically inertial and the black holes will move 
through it dynamically on near circular orbits. As a re- 
sult, the typical timescale for strong evolution of the met- 
ric at generic locations will be set by the orbital angular 
frequency 17, Torb ~ , and near each hole's horizon 
it will be much smaller than this. To follow these evolu- 
tions stably will require very short timesteps. At ^ Torb, 
which can be computationally very expensive. 

A promising alternative, proposed by Brady et. al., is 
to use co-rotating coordinates and a wisely chosen time 
coordinate, in which the metric components evolve (ev- 
erywhere, including near the horizons) on the gravita- 
tional radiation reaction timescale, t„ 3> Torb- In such 
coordinates, the timesteps can be much longer than in 
the conventional coordinates, thereby reducing the com- 
putation time by a large factor. 

Brady et. al. have proposed a specific choice for 



the lapse and shift functions a and (3j of numerical rel- 
ativity that, they argue, will force the coordinates to re- 
main co-rotating and the metric components to evolve 
on the long timescale Trr. This choice is one that min- 
imizes, at each time step, the integral of the (positive 
definite) square of the time derivative of the 3-metric, 
/ 7ij7fc;7*'^7-''d'^S, over the computation's spatial coor- 
dinate grid. This minimization leads to a set of minimal 
strain equations (MSE) |Q] for the lapse and shift: 
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(1) 
(2) 



where Di is the 3-covariant derivative. Substitution of (|l|) 
in (^ yields a second-order linear differential equation for 
/3fc: 



where the differential operator is 



(3) 



and 



(4) 

(5) 
(6) 



The purpose of this paper is to show that the system 
(^) is strongly elliptic and thus admits a unique analyti- 
cal solution for given Dirichlet boundary data, on a suffi- 
ciently small compact set. If our results can be extended 
to global existence and uniqueness, then the results will 
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confirm the MSE as an alternative tool for numerical so- 
lutions to the IBBH problem. 

The paper is organized as follows: Section 2 presents 
basic notions and results from regularity theory for el- 
liptic systems of partial differential equations. Section 

3 derives necessary and sufficient criteria for the strong 
ellipticity of the MSE and presents a step-by-step algo- 
rithm that can be used to test strong ellipticity for an 
arbitrary slicing of an arbitrary spacetime metric. In Sec. 

4 we use our algorithm to verify strong ellipticity for the 
optimal slicing of the Thorne test-bed metric (which 
is an approximation to that of a relativistic, inspiraling 
binary star), and we then argue that for any inspiraling 
binary the MSE are likely to be strongly elliptic. The na- 
ture of the boundary data is briefly addressed. Section 

5 concludes with a summary and comments on future 
numerical relativity schemes using the MSE. 

II. ELLIPTIC SYSTEMS OF PARTIAL 
DIFFERENTIAL EQUATIONS 

First, we introduce the notion of strong ellipticity, and 
then briefly discuss the integrability of strongly^ elliptic 
systems. 

Consider the following second-order differential equa- 
tion for a vector field defined on some n— dimensional 
manifold M: 

Ml[D,D]u^ = S'[D], (7) 

where M^[D,£'] includes the second-order derivative 
terms and all of the remaining terms are collected in the 
'source' vector functional 5'. The above equation is said 
to be strongly elliptic if, upon formal substitution at each 
point on M. of Dj by a vector field , the matrix thereby 
obtained, M^[^, has the following property: the eigen- 
values are all positive and bounded away from zero for 
all and all points in the manifold M (see e.g. Q). 

The Dirichlct problem for a vector field with analytical 
components u'"'(x) : Vl C M" ^ M can be formally stated 
as 

Xu*^ = , X e f^, (8) 
u^- - /^ X e 917, (9) 

where L is an elliptic differential operator and are 
analytical components of a vector field / defined on the 
boundary d^l. The Neumann problem consists in speci- 
fying the first derivative of normal to rather than 



* There are weaker notions of ellipticity that still guarantee 
the existence of a weak (defined in a distributional sense) so- 
lution, but they may not correspond to a unique analytical 
solution. Strong ellipticity avoids this ambiguity. 



itself. We shall be concerned with the Dirichlet prob- 
lem for the remainder of this paper, since our boundary 
data will consist solely of fit itself, (see, however, subsec. 
4. B). 

It can be shown (see e.g. [Q) that a strongly elliptic lin- 
ear system of second-order differential equations with an- 
alytical coefficients has a unique — defined on a compact 
subset W C n — analytical solution, for smooth Dirich- 
let boundary data. This establishes the well-posedness 
of the Dirichlet boundary value problem for the MSE, 
on sufficiently small regions, provided they are strongly 
elliptic. 

III. CONDITIONS FOR STRONG ELLIPTICITY 
OF THE MINIMAL STRAIN EQUATIONS 

A. Derivation of a Strong Ellipticity Condition 

Substituting Dj in M^'^ , as given by (Q), we ob- 

tain 

M> EE -x-'K^^K'^^eCi + \ (6^' + ^Hj^') ■ (10) 

When a basis has been chosen, this becomes a (3 x 3) real 
matrix. At any given point on the spacelike slice, we can 
introduce a local orthonormal basis, such that 7^ = 5ij. 
In this basis, M^^ becomes manifestly symmetric: 

(11) 

We can then rotate our basis, so as to diagonalize Kij-. 

Kij = dia.g{Kii,K22,K:i:i) = diag(fci, fc2, fcs). (12) 
Then Mik takes the form 

M,k = -X''hkkC^^k + \ + kSf) , (13) 

where there is no sum over i and fc, and the "==" denotes 
equality in the orthonormal basis in which Kij is diago- 
nalized. This can be further simplified by (i) introducing 
a unit-length vector field 

U = hx-^'\ (14) 

(where "unit-length" means ll + I2 + % — 1) and (ii) 
noting that Mij- is quadratic in ^i, which allows for the 
normalization 

lcT = E^.'-i, (15) 

i=l 

without loss of generality. Then Mik becomes 
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2M,fc 4 5,k + (1 - 2klk)C.^k- 



(16) 



B. Algorithm for Testing Strong Ellipticity 



The problem is now to determine what constraints on the 
vector field I (the normalized diagonal components of the 
extrinsic curvature in the orthonormal frame) ensure that 
the eigenvalues of Mik are positive and finite for arbitrary 
unit vector fields ^. 

By inspection, we see that one such constraint is that 
at most one of the U components vanish, as the following 
counter-example shows: if Z = ^ = (1, 0, 0), then 2Mik = 
diag(0, 1,1). It turns out that this is the necessary and 
sufficient condition, as we will now show. 

Let us consider an arbitrary vector field lu and contract 
it twice with Mik'. 



i,k 



i.k 
i,k 

>0, (17) 



where we have used the fact that = \l\ = 1. Since u) 
is arbitrary, this implies that the eigenvalues of Mik are 
non-negative. The inequality can be made strict if we 
assume that at least two of the li components are non- 
vanishing; then, since |/p — 1, we have \li\ < 1 for all i, 
whence 

u.\ < 161 ^ \hiMj^'^'\ < (18) 

and thus 

^M,fcw^w'=>0. (19) 

Therefore, the constraint that at most one of the k com- 
ponents vanish, is a sufficient condition for the strong 
ellipticity of the MSE. To show that this is also a nec- 
essary condition, we need to prove the converse: if the 
system (||) is strongly elliptic, then at most one of the li 
vanishes. Equivalently, we need to prove that, if two of 
the li vanish, then (^) is not strongly elliptic. 

Suppose, for definiteness, that h = h = 0; so h — 
±1. Then, choose the unit vector ^ to have components 
= 1, ^2 = 'Cs = 0- Equation ( p^ ) then implies that 
Mik is diagonal with entries Mu — 0, M22 — M33 = 1, 
which has a vanishing eigenvalue, and this implies that 
the MSE are not strongly elliptic. ■ 

Summarizing, the necessary and sufficient condition for 
the strong ellipticity of the MSE is that at most one di- 
agonal component of the extrinsic curvature — evaluated 
in an orthonormal frame in which Kij is diagonal — 
vanishes. 



The above strong ellipticity condition is only valid in 
an orthonormal basis in which the extrinsic curvature is 
diagonal. Testing this condition for some chosen met- 
ric and slicing requires a transformation to such a non- 
holonomic basis, which can be computationally tedious. 
It is therefore desirable to have a coordinate invariant al- 
gorithm for testing the strong ellipticity condition, which 
can easily be implemented for a general metric and ex- 
trinsic curvature in any basis one wishes. In this sub- 
section we derive such an algorithm. Our algorithm is 
expressed in terms of the object 



with trace 



JC = JCl^ —Kl 



(20) 



(21) 



In the orthonormal basis, 

K? ^ ih +I2 + hf = 1 + 2ihl2 + hh + hh). (22) 

Clearly, if /C^ 1 then at least two of the U components 
are non-zero, which implies that the MSE are strongly 
elliptic. This is a sufficient but not necessary condition, 
since the MSE can be strongly elliptic and JC^ = 1: the 
equation 



i=l 



= 0, 



(23) 



together with the constraint = 1' yields a 1- 

parameter family of solutions that geometrically corre- 
spond to two planes intersecting the unit sphere along 
two spherical triangles (which are a reflection of each 
other through the origin), whose sides are great circles. 
Except for the six corners of those triangles, where two 
li's vanish and therefore the MSE fail to be strongly el- 
liptic, along the edges the strong ellipticity condition is 
always satisfied. 

Further specialization can be accomplished by consid- 
ering the tensors 



KiklC) ± /Ci 



In the orthonormal basis, 

H± ^ diag(Z?, ID ± dia.g{hMM)- 



(24) 



(25) 



One of these vanishes if and only if two of the li compo- 
nents vanish, i.e., if and only if the MSE are not strongly 
elliptic (SE). 

Based on these results, we propose the following simple 
algorithm to test the strong ellipticity of the MSE in an 
arbitrary basis: 
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[1] X = ? Yes ^ I MSE are not SE | No ^ [2] 
[2]/C2 = 1 ? Yes ^ [3]. No 
[3] Is one of = ? Yes 



MSE are SE 



MSE are not SE 



No 



MSE are SE 



We note that an equivalent, coordinate invariant, crite- 
rion for strong ellipticity is simply the requirement that 
at most one of the eigenvalues of Kij vanishes, in any 
basis. This follows directly from the fact that a change 
of basis changes Kij into a similar^ matrix, thereby pre- 
serving the eigenvalues. 



IV. TESTING THE THORNE TEST-BED 
METRIC 

Thornc has constructed an analytic metric that he 
believes is similar, in its geometric and dynamical fea- 
tures, to those for inspiraling, fully relativistic neutron 
star binaries. He has suggested that this metric be used 
as a test bed for proposed lapse and shift functions for 
the IBBH problem. In this section, we ask whether the 
MSE equations for the lapse and shift are strongly elliptic 
in the Thorne test-bed metric. 

In co-rotating coordinates, {T,X,Y, Z}, the Thorne 
metric reads Q 

ds^ = -A^dT^ + TijidX' + B'dT){dX^ + B^dT), (26) 

where the co- rotating lapse A, shift , and 3-metric Tij , 
are exphcit functions of the spatial coordinates {X, Y, Z}, 
the two stars' initial mass m and their time-evolution sep- 
aration a(T); a(T) is freely adjustable but for approxi- 
mate radiation-reaction-driven inspiral can be chosen to 
be 



a{T) = flo 1 



1/4 



32m3r 
5ao 



The functions appearing in the metric ( |2q ) are: 

A=l- mF+, 
F± = [{X - af + Y'^ + Z^ + b^^ 

± [{X + af + Y^ + Z^ + b^] 
B = nY{P-l)dx 



1-1/2 



-1/2 



X -P X 



(27) 

(28) 
(29) 
(30) 



^Any two (n x n) matrices, A and B, are said to be similar 
if there exists a (n x n) matrix T, such that A — T^^BT. 
Similar matrices represent the same linear transformation, in 
different coordinate systems. 



p 


Ama 

~ (/?2 + a2)3/2' 




(31) 








(32) 


i 

y 


Sijil + mF+y, 




(33) 


TT 
ij 1 


^ h+{9 ® e - (j) ® (j)) + 


ft.x(^«)0 + 0«'^), 


(34) 


h+ 


= (1 + cos"^ t 

R 


1) cos[2<^- 217(r- i?)], 


(35) 




-Ama^VL^ , - 

= 2cos6'[20 

it 


-2f^(T-i?)], 


(36) 



R 



VR^ ~ Z^ 



R 



Ry/R^ - Z2 

1 



Vi?2 - ^2 i?2 _ ^2 

ndT ^(b + ^T. 



[Xdy - Ydx). 

{Xdv -Ydx), 



(37) 
(38) 
(39) 



Brady, Creighton and Thorne (BCT) have shown that 
for any metric (including Thome's test-bed) that has the 
form appropriate to a slowly inspiraling binary, there ex- 
ists a solution to the MSE that keeps the coordinates co- 
rotating to within an accuracy ~ TVot/Trr ~ i^Tn-)^^ <C 
1. If the MSE are strongly elliptic, and if appropriate 
boundary conditions are imposed on them, then that 
BCT solution should be their unique solution. The key 
issue of strong ellipticity should then be tested in the 
Thorne metric for a slicing that is close to that of the 
BCT solution. The natural such slicing is one of con- 
stant coordinate time t = T, since Thome's coordinates 
are co-rotating; and for that slicing the extrinsic curva- 
ture is 



^ 2A 



(40) 



T=To 



We have performed such a test, using our algorithm 
(end of Sec. 3), for the parameter set {M = 1.0, b = 
2.5, ao = 4.0}. Due to the symmetries of the problem, 
only one octant of R'^, locally defined by the cartesian 
basis {dx,dY,dz}, needed to be considered; we chose 
the region : < X,Y, Z < 10 and used a grid-size 
interval of 10""^. 

We found that 7^ for the points tested, but 
K,^ — 1 for a set of points topologically equivalent to a 
2-sphere with equator on the X — Y plane: for a given 
value of Z, the 2-surface given by IC^{X, Y) — l intersects 
the X — Y plane along a topological circle, whose radius 
decreases with increasing values of \Z\. For that set of 
points, Tifj ^ 0, and thus the MSE are strongly eUiptic. 
The cube extends sufficiently far from the two masses, 
that Kij and the other tensors built from it have the same 
behavior as at asymptotic spatial infinity; in particular, 
the relevant scalar quantities do not change sign. 
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A. Generality of the MSE's Strong Ellipticity 

We have shown that the MSE are strongly elliptic in a 
particular test-bed metric, which is an approximation to 
that of a relativistic inspiraling binary. It is reasonable 
to expect that for any realistic metric for inspiraling bi- 
naries, the MSE will also be strongly elliptic, since the 
failure of strong ellipticity requires a highly non-generic 
situation — at least two of the diagonal elements of Kij 
must vanish simultaneously in an orthonormal frame in 
which Kij is diagonal. Such non-genericity is likely to re- 
sult only from a very great spacetime and slicing symme- 
try and/or very special algebraic properties of the metric 
functions. 

A trivial example that illustrates the connection be- 
tween the failure of strong ellipticity of the MSE and a 
high degree of symmetry is the Schwarzschild metric: 



ds^ = -f{r)dt^ + r\r)dr^ + r^{d0^ + sin^ 
/(r) = 1 - 2Mr-\ 



(41) 
(42) 



With the usual Schwarzschild slicing (surface of constant 
time t), the MSE are not strongly elliptic: 



2V7(^ 



(43) 



In fact, for any static metric, dtgab — with slicing 
t = const., the extrinsic curvature vanishes identically, 
and the MSE are not strongly elliptic [cf. Eq. (^]. 
On the other hand, if we choose a Lemaitre slicing of 
Schwarzschild (corresponding to the reference frame of 
radially infalling observers): 



ds2 



-dt 



72 



2M 

r 



-df^ +r^{de^ 
di =dt + ^2Mr-^r^{r) dr. 



df = dt+ v/r/(2M)/-i(r) dr. 



r{r,i) 



/2M(f - i) 



2/3 



(44) 

(45) 
(46) 

(47) 



then the MSE are strongly elliptic. In particular, 



2M 



M 



,j - --dai]if,t) = y — diag(^,?-,rsin6l), (48) 

is positive definite everywhere except at the central sin- 
gularity r = 0, where Kff diverges and the other two di- 
agonal elements vanish. It is trivial to check that x 7^ 
and K.^ ^ 1, for all points (including r — 0, where both 
scalars diverge), thereby verifying strong ellipticity. 

In the Thorne test-bed metric, which has more sym- 
metries than a generic metric for inspiraling binaries, the 
condition JC^ ^ 1 (end of Sec. 3) failed to be obeyed 
only by a set of measure zero, for which the MSE were 
in fact strongly elliptic (the next necessary and sufficient 



condition for ellipticity was obeyed). This provides an 
example of how atypical a spacetime point would need 
to be in order for the MSE not to be strongly elliptic, 
even in an idealized metric. 

Based on the above observations, we conjecture that 
in a generic metric for inspiraling binaries, with a slicing 
near that which makes gab evolve on the timescale Trr, the 
MSE will be strongly elliptic, thence leading to a unique 
analytical solution for the choice of lapse and shift. 



B. A Note on the Nature of the Boundary Data 

On every spacelike slice, Et, appropriate boundary 
conditions for Pi need to be specified on an outer bound- 
ary 9Et that contains the two horizons. Hi and H2, in its 
interior, and on these two inner boundaries. Hi and i?2- 
Provided 5E( is sufficiently far from the two holes, we can 
take f3i on 9St to be the co-rotating one at asymptotic 
spatial infinity (e.g. like the one contained in Thome's 
test-bed metric ||^), which is a Dirichlet-type boundary 
condition. 

The inner boundaries are strong-field regions, and in 
general Neumann-type boundary conditions are needed 
1^,^ to ensure that /Si is indeed a minimal distortion one, 
which will typically require /?i ^ on Hi and H2 ■ Fixing 
the value of /3i itself, rather than its normal derivative 
on the inner boundaries, might lead to the pilling up of 
distortions near the boundary Unlike the Dirichlet 
problem, the Neumann problem does not admit a unique 
solution, even in an arbitrarily small compact set [cf. 
eqs. (8)-(9); we can always add to u'^ any function that 
is constant on each connected component of il], which 
precludes the a priori success of the minimal strain pre- 
scription for such type of boundary data. 

While this is still an open issue that deserves further 
investigation (a detailed understanding of how the met- 
ric evolves near the horizons through the IBBH phase), 
it may be that Dirichlet-type data may actually succeed 
in providing realistic and numerically stable boundary 
data on the inner boundaries. As long as adiabaticity 
is maintained, the shift generated by the MSE is such 
that the coordinates are kept co-rotating and the metric 
components evolving on the slow gravitational radiation 
reaction timescale, even near the horizons Q| . It is then 
reasonable to expect the normal ('time') derivative of the 
shift near the horizons to evolve rather slowly. In such 
case, we could then compute the inner boundary values 
of (3i on St from those on St_At (e.g. by changing them 
linearly), thereby avoiding having to solve a Neumann- 
type problem for the inner boundaries at every spacelike 
slice. For the initial slice, St=o, we could take f3i on the 
inner boundaries to be the co-rotating one, constructed 
from a post-Newtonian approximation for the early in- 
spiral phase. We must stress that this is a rather naive 
argument, and that only a detailed numerical study of 
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the spacetime evolution near the horizons, through the 
IBBH phase, can resolve this question unequivocally. 



V. CONCLUSIONS 

The BCT proposal for the optimal choice of lapse 
and shift in a 3 + 1 numerical relativity scheme is en- 
coded in a coupled set of second-order partial differential 
equations: the MSE. We have shown that the MSE are 
strongly elliptic if and only if at most one of the diago- 
nal components of the extrinsic curvature, evaluated in 
an orthonormal frame in which the extrinsic curvature 
is diagonal, vanishes; or, equivalently, if and only if the 
extrinsic curvature Kij possesses at most one zero eigen- 
value. An equivalent coordinate invariant algorithm was 
derived for general applicability to arbitrary choice of ba- 
sis and slicing. 

We tested this algorithm in the Thorne test-bed metric 
[H and found the MSE to be strongly elliptic. The fail- 
ure of strong ellipticity requires extremely atypical situa- 
tions, normally associated with highly symmetric metrics 
and/or special algebraic properties of the metric func- 
tions, which are very unlikely to occur in any realistic 
metric for inspiraling binaries. We therefore expect the 
MSE to be strongly elliptic for a generic metric for inspi- 
raling binaries. 

Evolving the IBBH spacetime in a 3 -I- 1 numerical 
scheme, using the minimal strain prescription, will typi- 
cally require solving a hyperbolic system for "fij and Kij , 
evolving forward in some time t, together with an ellip- 
tic system for the lapse a and shift (3i, at each surface 
of constant t, Ej. The main steps of such a computa- 
tion can be outlined as follows: (i) specify initial values 
for 7y and Kij] (ii) solve the elliptic MSE system (|^), 
using appropriate boundary values for Pi, to obtain (3i 
everywhere on that spacelike slice; (iii) solve Eq. (|l|) for 
the lapse a, using the j3i obtained from the MSE; (iv) 
solve the dynamical (hyperbolic) equations for ^ij and 
Kij using a and (3i produced by the MSE. This fully de- 
termines the spacetime geometry at the next spacelike 
slice, Et+At; (v) repeat steps (ii)-(iv), until the gravita- 
tional radiation reaction becomes too large for the evo- 
lution to be adiabatic. A simple, computationally in- 
expensive, self-consistency test would be to check that 
\\lij{t)/^ij[t — At) I — 1| ~ (JItit)^-^ ^ 1, for a given set 
of points at every spacelike slice. 
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